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I review recent work ]lj-|6j on magnetic dynamics of the high temperature superconductors using a 
model that combines two weakly interacting species of low-energy excitations: the antiferromagnetic 
spin waves which carry spin-1 and no charge, and Fermi- liquid- like quasiparticles which carry spin- 
1/2 and charge e. The model allows conversion of spin waves into electron-hole pairs; however, the 
low-energy spin waves are not collective modes of the quasiparticles near the Fermi surface, but rather 
are a separate branch of the low-energy spectrum. With certain experimentally justified assumptions, 
this theory is remarkably universal: the dependence on the detailed microscopic Hamiltonian and on 
doping can be absorbed into several experimentally measurable parameters. The z = 1 theory of the 
insulators and z = 2 theory of the overdoped materials, are both reproduced as limiting cases of the 
theory described here, which predicts that the underdoped materials remain in z = 1 universality 
class at sufficiently high temperature. This theory provides a framework for understanding both the 
experimental results and microscopic calculations, and in particular yields a possible explanation of 
the spin gap phenomenon. I also discuss some of the important unresolved issues. 

To appear in the Proceedings of the Stanford Conference on Spectroscopies in Novel Superconductors 
(1995). 



I. INTRODUCTION 



The evolution of the normal state properties of the high 
temperature superconductors with doping is schemati- 
cally described by the phase diagram in Fig.|l|. The 
stoichiometric insulators La2Cu04 and YBa2Cu306 un- 
dergo an antiferromagnetic transition at T/v, and exhibit 
short-range antiferromagnetic correlations well above 
Tjy. Upon doping by more than several per cent stron- 
tium or oxygen, these systems become metallic and no 
longer have a long-range antiferromagnetic order, but 
short-range magnetic correlations remain. It is univer- 
sally agreed that the normal state properties of this 
metallic state are far from the conventional metallic be- 
havior, and the term "strange metal" has been coined 
to describe them. Some of the key properties of the 
"strange metal" are: (i) short-range antiferromagnetic 
correlations, seen by NMR and neutron scattering; (ii) 
both the resistivity and the uniform magnetic susceptibil- 
ity increase as the temperature increases; (iii) some, but 
not all, of the materials exhibit a suppression of the low- 
frequency spectral weight for magnetic excitations for 
T < T* ~ 150-ftT (spin gap). At higher doping, the nor- 
mal state properties become increasingly similar to the 
Fermi liquid (FL) behavior with short-range antiferro- 
magnetic correlations. While Landau Fermi liquid theory 
in the orthodox sense requires temperature-independent 
spin correlations, the temperature dependence is weak in 
YBa2Cu307 and other fully doped materials, where close 
proximity of the FL state is evident. 
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FIG. 1. A phase diagram of the high-temperature su- 
perconductors as a function of doping and temperature. In 
the limit of zero doping, the insulating parent materials de- 
velop a long-range Neel order and their low-energy excita- 
tion spectrum has only one species of excitations — the spin 
waves. The overdoped materials at low temperatures exhibit 
Fermi-liquid-like properties; there, the quasiparticles are the 
only species of low-energy excitations. We attempt to un- 
derstand the properties of the intermediate "strange metal" 
phase, where superconductivity occurs, by combining the ex- 
citations found in the extreme doping limits, namely the spin 
waves and the quasiparticles, and allowing spin waves to con- 
vert into electron-hole pairs. 

High-temperature superconductivity occurs in the in- 
termediate strange metal phase in Figj^. Currently, there 
is no consensus regarding the precise microscopic Hamil- 
tonian of this phase. In what follows, I describe a theoret- 
ical approach to this problem which bypasses the missing 
information about the microscopic Hamiltonian by as- 
suming continuous evolution of the low-energy magnetic 
and charge excitation spectrum between the antiferro- 
magnetic insulator and the Fermi liquid phases, through 
the intermediate strange metal phase (see, e.g. and 
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references therein, and relevant earlier work M [L(J ) . 

II. NEARLY ANTIFERROMAGNETIC FERMI 
LIQUID 

The approach of Refs. is based on combining the 
low-energy excitations encountered in the two extreme 
doping limits in Fig.jjJ the undoped antiferromagnetic in- 
sulators where the low-energy excitations are spin waves, 
and the overdoped limit where the low-energy excitations 
are the usual Fermi liquid quasiparticles. I emphasize 
that such an approach makes no explicit assumptions 
about the microscopic Hamiltonian, and also it makes 
no assumptions regarding the mechanism through which 
the two types of low energy excitations are formed. In- 
stead, the theory is based on the premise that at low 
energies, the spectral weight of the system is shared be- 
tween the spin waves and the quasiparticles, and that the 
primary interaction mechanism is spin wave to electron- 
hole pair conversion. A microscopic basis for this theory 
is provided by the mean-held spin density wave (SDW) 
calculations, described in ||. 

The presence of the first species of low-energy excita- 
tions, the quasiparticles, is evident from photoemission 
experiments. Such quasiparticles at low energies must be 
similar to the quasiparticles of the Landau Fermi liquid 
theory, and carry both charge e and spin S = 1/2, in or- 
der to account for the sharp quasiparticle peak near the 
Fermi energy, observed in photoemission. The quasipar- 
ticles are gapless and form a Fermi surface, the shape of 
which is very important for the interaction between the 
quasiparticles and spin waves. 

The other species of low energy excitations, the spin 
waves, are the only low-energy modes in the antiferro- 
magnetic insulators La2CuC>4 and YBa 2 Cu30g. Below 
Tyv, the spin waves are gapless, and in small-g limit have 
a linear spectrum cu q ~ cq and a mean-free path that is 
longer than the wavelength. When elevated temperature 
or doping destroy the long-range order, and only short- 
range antiferromagnetic correlations remain (the corre- 
lation length £ is finite), the spin waves are no longer 
gapless and also are overdamped for q < Never- 
theless, for larger wavevectors q > the spin waves 
remain nearly the same as in the presence of long-range 
order, because they sample only the local order for dis- 
tances of order wavelength A = 2ir/q <C £, and therefore 
are nearly insensitive to the absence of antiferromagnetic 
correlations at distances larger than £. The evidence for 
their existence in the strange metal phase is based pri- 
marily on NMR and neutron scattering measurements, 
and is discussed in what follows. The spin waves carry 
spin 5=1 and no charge. 

Whereas on the microscopic level an antiferromagnetic 
spin wave can be regarded as a coherent electron-hole 
pair, the contribution from the quasiparticle states near 



the Fermi surface accounts for only a fraction of the to- 
tal spectral weight of the spin wave, if there is no nest- 
ing. The spin waves are therefore a separate branch of 
excitations, rather than a collective mode of low-energy 
quasiparticles. In that sense, the theory described here 
is somewhat similar to the spin-charge separation theory, 
where low-energy spin excitations cannot be represented 
as collective modes of low-energy charge excitations, even 
though on the microscopic level both originate from elec- 
trons. The difference between our model and the 
spin-charge separation picture applied to the high tem- 
perature superconductors Jul, is that in our model one 
species of excitations (spin waves) carries only spin, while 
the other (quasiparticles) carries both spin and charge, 
therefore there is charge separation, but not spin sep- 
aration. Recently, Sachdev (|] and Laug hlin H have 
pointed out that in a system with spin-charge separa- 
tion, some of the spinon and holon excitations may form 
bound states which carry both spin and charge, and in 
many respects are similar to conventional quasiparticles. 




FIG. 2. Magnetic excitation spectrum of type (B) model 
(classification due to Millis Jl^]), where Q = (n/a,n/a) can 
connect two points on the Fermi surface. The minimum of 
the spin wave dispersion u) q = (c 2 |q — Q| 2 + A 2 ) 1//2 is located 
inside the electron-hole continuum. 

One of the key factors that affect the magnetic and 
charge properties is the position of the minimum of the 
spin wave spectrum, Q = (w/a,iv/a) for the commensu- 
rate short range order discussed here, with respect to the 
electron-hole continuum. If the spin wave spectrum is 
located inside the electron-hole continuum, as it is illus- 
trated in Fig.^, the spin waves can convert into electron- 
hole pairs. This conversion process becomes the domi- 
nant contribution to spin wave damping at low energies, 
and has a profound effect on the magnetic properties. 
In the classification introduced by Millis, this situation 
is called model B. For a discussion of the opposite case, 
when Q cannot connect two points on the Fermi surface 
(model A), see Refs. [0|J|. 
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The results of recent photoemission measurements |13|] 
indicate that model B is relevant to the Y- and Bi-based 
high-temperature superconductors. In what follows, I 
concentrate exclusively on model B. 

III. SPIN GAP 

The spin gap behavior is the suppression of the low 
frequency spectral weight for magnetic excitations, which 
can be observed by NMR and neutron scattering. In the 
model at hand, this suppression is expected already in 
a very crude approximation which treats spin wave and 
quasiparticle excitations as noninteracting, free modes. 
Below, I describe a qualitative picture of the spin gap 
behavior suggested by Sokol and Pines |§Q , and derived 
starting from model B by Sachdev, Chubukov, and Sokol 

In the non-interacting approximation, the density of 
states for magnetic excitations is a sum of spin wave 
and quasiparticle contributions. For frequencies below 
the gap for spin wave excitations, u> < A = c/£, only 
the quasiparticles contribute to the magnetic density of 
states. For frequencies to > A, both the spin waves 
and the quasiparticles contribute. In a measurement of 
the total magnetic density of states, for example in a 
NMR relaxation or neutron scattering measurement at 
low temperatures, one expects to see a gap-like suppres- 
sion of the density of states for magnetic excitations for 
w < A. 

Beyond the non-interacting approximation, the life- 
time of spin waves primarily is due to their conversion 
into electron-hole pairs. Upon increasing doping, the spin 
gap is washed out by a combination of two separate ef- 
fects: first, the spin waves become overdamped because 
the rate of conversion increases; second, the gap for spin- 
wave excitations A = c/£ increases as £ decreases. The 
remaining spin waves, which are pushed to much larger 
energies, can in principle be seen by neutron scattering. 
At least two recent experiments seem to allow such an in- 
terpretation: the evidence for the magnetic character of 
the 41 meV peak in YBa2Cu3C>7 by Keimer et al. flif] 
and the observation of the zone-boundary magnon in 
La2- x Sr x Cu04 by Aeppli et al. jl5|]. Keimer [|14j and 
Barzykin and Pines || have suggested that the 41meV 
feature is a coherent, triplet, electron-hole pair - in other 
words, a spin wave (several alternative interpretations 
have been suggested as well). The zone boundary high- 
energy spin wave (a magnon), which continuously evolved 
from the identical mode in the insulator, seems to be 
most likely explanation of Aeppli's measurements; how- 
ever, this data is very preliminary Jig]. 



IV. UNIVERSAL THEORY OF MAGNETIC 
DYNAMICS 

In this section, I describe recent work by Sachdev, 
Chubukov, and Sokol ||, where the universal behavior 
of the magnetic correlations near Q, and of the bulk sus- 
ceptibility, was calculated using the rate of the Landau 
damping of spin waves as an input parameter. Several 
important assumptions about the microscopic model are 
built into this theory. First, we limit our study to sys- 
tems with short-range order at T — 0. The region at low 
doping where the long range Neel order exists at T = 
requires separate consideration. Second, we require that 
there is no nesting, i.e. the areas of the Fermi surface that 
are adjacent to the points connected by Q, are not paral- 
lel to each other. The photoemission measurements |u| 
show that this is the case in Y- and Bi-based supercon- 
ductors. Third, we assume that there is no substantial 
q-, lo— 1 and T— dependence of the Landau damping for 
q — Q ~ l/£ and u> ~ A, respectively. This assump- 
tion is more difficult to verify experimentally; it holds in 
mean field theory on the disordered side if the correlation 
length is large enough. 

Under these assumptions, the additional damping due 
to the spin wave to electron-hole pair conversion can be 
included into the theory by inserting its rate T into the 
noninteracting spin wave dynamical response function (a 
similar expression was introduced by Barzykin, Pines, 
Sokol, and Thelen |[6[ on phenomenological grounds): 

const const 
X q,w -> — — , (1) 

where the spin wave spectrum for finite £ (short-range 
correlations) has the following form: 

uj q = Vc 2 |q-Q| 2 + A 2 , A = c/t (2) 

The exact dynamical susceptibility of the model is ob- 
tained by calculating the effects of mutual scattering of 
spin waves, primarily the additional damping, on x(q, w ) 
given by Eq.(|l|). When the Landau damping processes 
dominate dissipation, the result of such a calculation dif- 
fers only slightly from Eq. ([[]). 

With these assumptions, the dynamical magnetic sus- 
ceptibility near the staggered wavevector, q Q <C a -1 , 
is given by: 

x(q>w) = — {-) $ s I — - — ,-J, (3) 

and the temperature-dependent part of the uniform mag- 
netic susceptibility, Xm by: 

Xu(T) - X u(T = 0) = (1 + a')^ u . (4) 

Here, a' is a non-universal constant, both Q s and $ u are 
universal and computable functions of their arguments, 
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and T] is a universal critical exponent which is very small 
and for all practical purposes can be replaced by zero. 

Eqs.(^)J^) apply when all dimensionful parameters and 
variables are smaller than the respective lattice cutoffs, 
which roughly translates into the following: 



w,T,A,T<^min(J,E F ), f -1 , |q- Q| « a" 



(5) 



where J is the exchange constant, Ep the Fermi energy, 
and a the lattice spacing. If the conditions (|J) apply, 
the magnetic dynamics near Q depends on four dimen- 
sionful parameters, two of which describe the spin wave 
spectrum (c and A), one is set by the size of spin (Z), and 
one is determined by the rate of the Landau damping (r). 
The temperature-dependent part of the bulk susceptibil- 
ity, which adds to the Pauli term, is also universal upto 
the overall prefactor 1 + a'. 

For a low-energy theory this number of parameters is 
not excessive because they include all the information 
about the material that is reflected in its magnetic be- 
havior. Three of them are the overall scales of magnetic 
moment, distance, and energy, and are therefore unavoid- 
able in any low-energy theory; determination of these 
three parameters from the experimental data is straight- 
forward and is a matter of simple normalization. The 
fourth parameter serves as the single crossover variable 
which describes the evolution from the scaling behavior 
of the insulators (the dynamical exponent z = 1) to that 
of the overdoped materials (z = 2). One of the key results 
of H is that there is only one such crossover parameter. 
Here the dynamical exponent z relates temperature de- 
pendences of the characteristic frequency scale for mag- 
netic correlations, Co, and the characteristic length scale 
(the correlation length, £): 
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FIG. 3. Evolution of x£(q,id) = / dqx"(q,w) at T = 
as a function of the dimensionless crossover parameter r/A. 
The underdoped materials YBa2Cu30e.63 and YE^CruOs 
are in regime (b), and have 2 = 1. The same applies to 
Lai.ssSro.isCuCU at high temperatures; at low temperatures, 
some modifications of the theory due to the incommensura- 
bility of short-range order are necessary. The fully doped 
material YBa2Cu307 is in regime (c), and has 2 = 2. 



Fig.|| illustrates the evolution of zero temperature mag- 
netic properties with doping on the example of the local 
dynamical magnetic susceptibility: 



J dqx"(cbw), 



(7) 



which shows a gradual crossover, observed as a function 
of r/A, from z = 1 regime (r/A < 1 §§,§) to z = 2 
regime (r/A > 1 P,pTj|). The crossover between z = 1 
and z — 2 regimes at low temperature is expected, and 
observed experimentally, as doping increases; in the uni- 
versal theory, a similar crossover is also expected at a 
fixed doping as the temperature increases; it is however 
not observed experimentally because scaling fails alto- 
gether when at high temperatures the correlation length 
becomes very short, £/a ;$ 2 (Barzykin and Pines [j6|). 
The z = 2 phase shows a further crossover to the Lan- 
dau Fermi liquid when A 2 /r becomes larger than u> or 
T. FigJ^ describes the sequence of crossovers observed as 
r/A increases. 



A. z=l limit (r/A < 1) 

The limit Y / A = and A > corresponds to an insu- 
lator with short-range antiferromagnetic correlations at 
T = (Fig.||a). Such an insulator has a true energy 
gap A because creating any combination of excitations 
(spin waves which have a gap) above the ground state 
requires a finite energy, hence x'Li 10 < A) = 0. While 
this limit has no direct relevance to the cuprate oxide in- 
sulators (which develop Neel long-range order at T = 
and have gapless spin waves), it helps to understand the 
behavior of the underdoped materials where r/A < 1. 
The effect of small doping is shown in Fig^b: x'L( w ) be- 
comes finite for all lo, but it remains much smaller for 
lu < A compared to to > A. The gap which existed for 
r/A = transforms into a knee-like feature at u> ~ A 
for r/A < 1. This behavior reproduces the essential fea- 
tures of the spin gap observed in some of the underdoped 
high-temperature superconductors. 

At high temperatures T ^> A, the behavior of r/A < 1 
underdoped system is similar to that of the Yj A = in- 
sulator, because for energies u> > A the spectral weight 
of the two systems is nearly the same, and in both cases 
is dominated by the spin wave contribution (compare 
Figs.||a,b). As a result, the magnetic dynamics in this 
temperature range is in the z = 1 universality class, and 
exhibits quantum critical behavior Co ~ A ~ T . 

The NMR measurements in the underdoped materials 
YBa 2 Cu 3 6 . 6 3 (Takigawa 0) and YBa 2 Cu 4 8 (Imai 
et al. |||, Stern et al. |TJ, and Corey et al. @) are 
consistent with z — 1 predictions by Sokol and Pines 
0]: both materials have a spin gap at low tempera- 
tures, and X1T/T2G ~ const at high temperatures. In 
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Lai.85Sr .i5CuO4, I/T2G nas n0 * been measured at 
sufficiently high temperatures (T — 500 — lOOOif), but 
1/Ti has been, and it is nearly the same as in the insula- 
tor (Imai et al. pl| ). Therefore, it has to be dominated 
by the spin wave contribution, which obeys the z = 1 the- 
ory. At low temperatures, the application of this theory 
to Lai.85Sro.i5Cu04 requires modifications due to the in- 
commensurability of short range order. According to the 
theory described here, Lai.gsSro.isCuC^ should have a 
smaller spin gap than the underdoped YBCO materials, 
which appears to be almost obscured by the supercon- 
ducting phase (see Ref. for a discussion). 



B. z=2 limit (T/A > 1) 

Figs.||c,d describe different regions of the same fre- 
quency dependence of at T — in the overdoped 
case, where T/A > 1: (c) corresponds to ui ~ A 2 /r, and 
(d) tow< A 2 /r. Regime (c) is z — 2 quantum critical; 
regime (d) is z = 2 quantum disordered, which is equiva- 
lent to the Landau Fermi liquid. The theory of magnetic 
behavior in the z = 2 regime has been developed by Mil- 
lis, Monien, and Pines || and Millis [{U^. 

For temperatures T ^ A 2 /r, the characteristic fre- 
quency scale for magnetic correlations Co ~ T modulo 
log(A/T) corrections, and the correlation length is given 
by Eq.@ with z = 2. In the opposite limit T <C A 2 /r, 
the correlation length is nearly T- independent, as is ex- 
pected for the Landau Fermi liquid. 

The fully doped material YBa2Cu3C>7 is in the z = 2 
universality class, and its behavior between T c and room 
temperature is intermediate between the z = 2 quan- 
tum critical regime (Figj^c) and z — 2 quantum dis- 
ordered regime (Fig.gjd). The NMR experiments find 
T1T/T2Q ~ const in agreement with the z — 2 quantum 
critical prediction (Imai et al. p2| ); however, both of the 
rates do not vary strongly with temperature, which is 
an indication of the proximity of the Fermi liquid phase, 
where they saturate. This material does not have a spin 
gap, in agreement with its assignment as being in the 
z — 2 regime. 



V. UNIFORM SUSCEPTIBILITY 

The approach of Ref. || allows one to calculate the 
temperature-dependent part of the uniform spin suscep- 
tibility, Xu(T) — Xu{T = 0), and correctly predicts that 
it monotonically increases as a function of T above T c , 
that its slope becomes smaller as doping increases, and 
that temperature-independent Pauli Xu(T) is recovered 
in the overdoped case. 
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FIG. 4. A sketch of the uniform spin susceptibility Xu(T) 
in the underdoped YBCO (solid line), and the theoretical 
prediction of Eq.(^) (dashed line). Since Eq.(^J) is only valid 
in the normal state, comparisons below T c are not meaning- 
ful. Above T c , the experimentally observed Xu{T) increases 
rapidly, and its slope is qualitatively consistent with that pre- 
dicted by Eq.(^j). At higher temperatures, the theory pre- 
dicts continuing linear increase, whereas the experimental 
data shows saturation at T ~ 150 — 200A". Following the 
analogy with similar behavior of Xu{T) for the insulating an- 
tiferromagnets at T ~ 0.7 — 0.9J, one may speculate that the 
discrepancy is due to lattice corrections. 

However, the experimentally observed slope of Xu(T) 
sharply decreases above T ~ 150 — 200K compared to 
the slope below this temperature but above T c |2^,[l7],|l9| , 
whereas the universal theory predicts that the increase 
of Xu{T) continues with the same slope (Fig.^). One 
may speculate that this deviation is due to lattice cor- 
rections; detailed calculations must be performed to de- 
termine whether or not this speculation is correct. 

Millis, Ioffe, and Monien [^4| recently presented a sce- 
nario of spin gap formation based on spin-charge sepa- 
ration, accompanied by singlet pairing of spin-carrying 
excitations (spinons). Since holons do not have spin and 
therefore do not contribute to the magnetic density of 
states, the bulk susceptibility must vanish upon singlet 
pairing of spinons at low temperatures. They argued 
that the observed decrease of the bulk susceptibility at 
low temperatures in underdoped YBCO can only be ex- 
plained with this scenario. 

However, their conclusion relies on the assumption that 
both x'Li 1 ^ ~* Oj^ 1 ) and Xo(T) would extrapolate to zero 
at low temperatures even without the superconducting 
transition, a statement that cannot be verified experi- 
mentally. In our view no conclusion can be drawn 
from the experimental data on Xu(T) above T c as to 
whether at T — > the magnetic density of states with- 
out superconductivity would decrease all the way down 
to zero, as it is required by their model, or would become 
much smaller than at high temperatures, but remain fi- 
nite, as it is expected in our model. 
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VI. UNRESOLVED ISSUES AND CHALLENGES 

The universal scaling theory yields a qualitative sce- 
nario for the evolution of magnetic behavior with doping 
which is consistent with the experiment. It also has been 
very successful in quantitatively explaining as well as pre- 
dicting the results of some of the experimental measure- 
ments, notably the NMR relaxation rates in the insula- 
tors and doped materials. However, it fails to explain 
the detailed temperature dependence of Xu(T) at high 
temperatures. This failure is likely to be caused by the 
lattice corrections in the broad sense: the influence of 
the Brillouin zone boundary on q-intcgrations, the non- 
linearity of the spin-wave spectrum at large wavevectors, 
bilayer coupling in YBCO, or energy dependence of the 
Landau damping rate T. 

Different quantities are differently affected by the lat- 
tice corrections, and some are more robust than others. 
For instance, over a range of temperatures, \u in the 
Heisenberg model is affected by the lattice corrections, 
while the characteristic frequency to is not. Furthermore, 
some of the lattice corrections that do appear, can be 
absorbed into the dimensionful parameters of the model, 
while the universal scaling functions remain nearly un- 
affected [§5|,|6| . A detailed study of these lattice effects 
may allow their inclusion into a phenomenological exten- 
sion of the scaling theory to expand its range of applica- 
bility. 

Another important challenge is to determine the pre- 
cise role of bilayer coupling, which is central to the sce- 
nario of spin gap formation by Millis and Monien p7[ . 
In the theory described here, the bilayer coupling affects 
the dimensionful parameters of the model, but does not 
affect any observables expressed as a function of these 
parameters, unless the size of bilayer coupling is compa- 
rable to other low-energy scales. Note that in our model, 
the spin gap may form with or without bilayer coupling, 
but it is enhanced if bilayer coupling is present. 

Finally, it is very important to understand on the mi- 
croscopic level how the spin waves and the quasiparticles 
share the spectral weight for intermediate doping, and in 
particular the mechanism of the experimentally observed 
reduction of the Pauli contribution to the uniform spin 
susceptibility in the underdoped case. The experimental 
data indicates that the spin wave to electron-hole pair 
conversion rate is also suppressed in the underdoped ma- 
terials, which in the theory described here results in spin 
gap formation. The cause of such a reduction is not fully 
understood, and may be a consequence of a reduced den- 
sity of states for the quasiparticles at low doping, which 
is also seen in Xu{T — 0). 

Studies of microscopic models for the high-Tc cuprates 
should help to clarify whether the experimentally ob- 
served evolution of the key parameters of the model dis- 
cussed here is similar to that obtained in microscopic 



calculations, thereby allowing one to check the validity 
of the proposed theory. 
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